Zitterbewegung of electronic wave packets in semiconductor nanostructures 
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We study the zitterbewegung of electronic wave packets in III-V zinc-blende semiconductor quantum 
wells due to spin-orbit coupling. Our results suggest a direct experimental proof of this fundamental 
effect, confirming a long-standing theoretical prediction. For electron motion in a harmonic quantum 
wire, we numerically and analytically find a resonance condition maximizing the zitterbewegung. 



The emerging field of spintronics has generated a large 
deal of interest in the effects of spin-orbit coupling in 
semiconductor nanostructures and their possible appli- 
cations [1]. Spin-orbit coupling is a relativistic effect 
described by the Dirac equation and its nonrelativistic 
expansion [2]. In semiconductors spin-orbit coupling of 
itinerant electrons is much stronger than in vacuum. This 
is due to the interplay of crystal symmetry and the strong 
electric fields of the atom cores [3]. In fact, the effec- 
tive models describing the band structure of III-V semi- 
conductors have many similarities to the Dirac equation. 
However, the fundamental gap between conduction and 
valence band in a semicondcutor is of order an electron 
volt or less, which is much smaller than the gap between 
solutions of positive and negative energy of the free Dirac 
equation. This observation can serve as a general heuris- 
tic explanation for the importance of spin-orbit effects in 
semiconductors [3]. 

Another prediction of relativistic quantum mechanics 
is the zitterbewegung of electrons [2], which, however, 
has not been observed experimentally yet. For free elec- 
trons, i.e. in the absence of an external potential, such 
an oscillatory motion occurs if solutions of both positive 
and negative energy of the free Dirac equation have a fi- 
nite weight in a given quantum state. In this letter we 
investigate the zitterbewegung of electron wave packets 
under the influence of the enhanced spin-orbit coupling in 
III-V zinc-blende semiconductor quantum wells. Strong 
spin-orbit coupling generally requires large gradients of 
the external potential, as they are provided by the heavy 
atom cores in such systems. 

An important effective contribution to spin-orbit cou- 
pling in such systems is the Rashba term which is due to 
structure-inversion asymmetry of the confining potential 
and takes the following form [4] 



H R = (a/h)(p x ay -p y a x ). 



(1) 



where p is the momentum of the electron confined in 
a two-dimensional geometry, and a the vector of Pauli 
matrices. The Rashba coefficient a is essentially pro- 
portional to the potential gradient across the well and 
therefore tunable by an external gate. Thus, the single- 
particle Hamiltonian is given by H = p 2 /2m + li.R, where 



m is the effective band mass. The components of the 
time-dependent position operator 



fff(t) = e lH/h r(0)e- lH/h 
in the Heisenberg picture read explicitly 
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where the operators p and a on the right hand sides are in 
the Schrodinger picture and therefore time-independent. 

We now proceed by evaluating the above time- 
dependent position operators within a Gaussian wave 
packet with initial spin polarization along the z-direction 
perpendicular to the quantum well, 
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Clearly we have (i^\r\tp) = 0, (if}\p\ijj) = hko, and the 
variances of the position and momentum operators are 
(Ax) 2 = (Ay) 2 = d 2 /2, (A Px ) 2 = (A Py ) 2 = h 2 /2d 2 . 
Thus, the group velocity of the wave packet is given by 
hko/m , while its spatial width is described by the pa- 
rameter d with the minimum uncertainty product typical 
for Gaussian wave packets, Ap x Ax — Ap y Ay — Ti/2. 
A direct calculation gives 

/.2tt 



(1>\x H (t)\1>) = 



O.r 



-d k 'i 



2tt 



dip sin tp 



o 



dq e 



-q 2 -\-2dq(k() x cos <p-\-ko y sin tp) 



(2aq 



(6) 



1 



In the above expression, ip is a usual polar angle in the 
xy-plane, and q is a dimensionless integration variable. 
The time dependence in the integral can be viewed as a 
zitterbewegung the electron performs under the influence 
of spin-orbit coupling. Clearly, this integral contribution 
vanishes for k 0y = 0, i.e. if the group velocity is along 
the ^-direction. More generally, one finds that 



(tp\k ■ f H (t)/ko\ip) = hk t/m. 



(7) 



which means that the zitterbewegung is always perpen- 
dicular to the group velocity of the wave packet. Let us 
therefore concentrate on the case ko x = 0. By expanding 
the exponential containing the trigonometric functions 
one derives 

(iP\x H (m = ^ r (l-e- d2k2 °y) 
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Thus, the amplitude of the zitterbewegung is propor- 
tional to the wave length of the electron motion perpen- 
dicular to it, and the oscillatory zitterbewegung changes 
its sign if the translational motion is reversed. If the 
product dkoy is not too large, dko y < 1, only low values of 
the summation index n lead to substantial contributions, 
and the Gaussian factor in the integrand suppresses con- 
tributions from large values of q. Thus, a typical scale 
of this integration variable is leading to sizable contri- 
butions is q w l/y/2. Thus, a typical time scale in the 
integrand is T = \/2irhd/a, and when averaging the zit- 
terbewegung over times scales significantly larger than T, 
the cosine term drops giving 



(V|zff(*)|V> = (1/2***) (1 - cM-d 2 kl y )) , (9) 

i.e. the time-averaged guiding center of the wave packet 
is shifted perpendicular to its direction of motion. Note 
that the zitterbewegung is absent for ko y — [5]. 

In the opposite case dko y 3> 1 the Gaussian approaches 
a (5-function. In this limit one finds (for ko x = 0) 

(i>\x H (tM) = (l/2k 0y ) (1 - cos (2ak Qy t/h)) . (10) 

Here the frequency of the zitterbewegung is Q = 
2ako y /fi, and the guiding center of the wave packet is 
also shifted in the direction perpendicular to its group 
velocity. Note that fifl is the excitation energy between 
the two branches of the Rashba Hamiltonian H at a given 



momentum k 
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The zitterbewegung of an electron in a quantum well 
as described above is naturally accompanied by a broad- 
ening of the wave packet, where the dominant contri- 
bution stems from the dispersive effective-mass term 



in the Hamiltonian. Such a broadening might pose 
an obstacle for experimentally detecting the zitterbewe- 
gung. However, the broadening can be efficiently sup- 
pressed and limited if the electron moves along a quan- 
tum wire. In fact, the motion of electrons in quantum 
wells is generally under better control if additional lat- 
eral confinement is present. We therefore consider a har- 
monic quantum wire along the y-direction described by 
H = p 2 /2m + mu) 2 x 2 /2 + Hr, where the frequency w 
parametrizes the confining potential perpendicular to the 
wire [6,7]. For this case exact analytical progress as above 
does not seem to be possible, and we therefore follow a 
numerical approach combined with an approximate ana- 
lytical study. To be specific, we consider an electron with 
a given momentum ko y along the wire and injected ini- 
tially into the lowest subband of the confining potential 
with the spin pointing upwards along the z-direction, i.e. 
the initial wave function for the a;-direction is a Gaussian 
whose width is determined by the characteristic length 
A = yjh/muj of the harmonic confinement. In our numer- 
ical simulations of the exact time evolution we find again 
a zitterbewegung perpendicular to the electron motion 
along the wire with the width of the wave function across 
the wire being limited by the confining potential. More- 
over, the amplitude of the zitterbewegung becomes max- 
imal if the resonance condition \KCl\ = \2ako y \ — frco is 
fulfilled. This general finding is illustrated in Fig. 1 where 
the wave number along the wire is fixed to be ko y X = 5 
and the Rashba parameter a is varied around the reso- 
nance condition. Equivalent observation are made if the 
Rashba coupling is fixed while the wave number ko y is 
varied. In Fig. 2 we have plotted the amplitude of the 
zitterbewegung as a function of tt/cu — 2ako y /hu> for 
different values of the wave number ko y along the wire. 
In this range of parameters, the resonance bocomes nar- 
rower with increasing fco y , while its maximum value is 
rather independent of this quantity and remarkably well 
described by \j\[2. 

A qualitative explanation for this resonance can be 
given by writing the Hamiltonian in the form H = 
Ho+Hi with Ho = huj(a+a+l/2) + h 2 k l y /2m + ako y a x , 

Hi = —iy/hmu)/2(a/h)(a — a + )a y , and a, a + being the 
usual harmonic climbing operators [7]. The zitterbewe- 
gung is induced by the perturbation Hi which can act 
most efficiently if the unperturbed energy levels of Ho 
are degenerate having opposite spins. This is the case at 
|2afco y | = fiuj. 

Another way to understand this resonance condition 
is to consider a truncated model where the Hamiltonian 
has been projected onto the lowest to orbital subbands, 
an approximation which is known to give very reasonable 
results for the low- lying energy spectrum of the wire [7] . 
For a given wave number ko y the truncated Hilbert space 
is spanned by the states |0, |), |0, J.), |1,T)j |1,|), where 
the arrows denote the spin state with respect to the z- 
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direction, and and 1 stand for the ground state and the 
first excited state of the harmonic potential, respectively. 
When applying the transformation 
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the projected Hamiltonian and in turn the time evolution 
operator become block-diagonal, 



Ue -tm u+ = (h + {t) \ 
I h-(t)J 



(12) 
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and fl± — (±a^Jmuo /2h, 0, — Huj/2 =F ak ny ). Now using 

one obtains for the initial state |0, f) the following time- 
dependent expectation value 



(o,t|^(i)|o,t) = 



V2 M+M- 

• sin (fi + t/~h) sin (/x_i/?i) . (15) 



At resonace, Ulo = ±2ako y , we have /i^ = 0, and if 
can be neglected compared to /zfj_ (which is the case for 
large enough ko y ) the amplitude of the zitterbewegung is 
approximately X/V2, in accordance with the numerics. 

We note that the above resonance condition can be 
reached by tuning either the Rashba coupling or the wave 
vector of the electron injected into the wire. For instance, 
for a wire of width A = 50nm and the band mass of InAs 
being 0.023 in units of the bare electron mass, one has 
fiui s=a lmeV. Typical values for the Rashba parameter a 
in InAs are of order 10~ n eVm [8-14] leading to a wave 
length Aoj, = 2ir/ko y ~ lOOnm at resonance. For GaAs, 
the band mass is larger while the Rashba coefficient is 
typically an order of magnitude smaller than in InAs [15], 
giving a wave length of a few ten nm at resonance for 
A = 50nm. 

We are propose that electron zitterbewegung in semi- 
conductor nanostructures as described above can be ex- 
perimentally observed using high-resolution scanning- 
probe microscopy imaging techniques as developed and 
discussed in Refs. [16,17]. As a possible setup, a tip suffi- 
ciently smaller than the width of the wire can be moved 
along the wire and centered in its middle. For an ap- 
propriate biasing of the tip, the electron density at its 



location is depleted leading to a reduced conductance of 
the wire. Since the amplitude of the zitterbewegung re- 
flects the electron density in the center of the wire, the 
zitterbewegung will induce beatings in the wire conduc- 
tance as a function of the tip position. These beatings are 
most pronounced at the resonance, see Fig. 1. Note that 
the oscillations shown there as a function of time can be 
easily converted to the real-space y-coordinate by mul- 
tiplying the abscissa by Hkoy/m. Generally we expect 
spin-orbit effects in STM experiments to be more pro- 
nounced in the presence of additional confinement such 
as in a quantum wire. 

Moreover let us consider spin-orbit coupling of the 
Drcsselhaus type [18-20], 



(16) 



which is due to the bulk-inversion asymmetry in zinc- 
blende semiconductors, and the coefficient [3 is deter- 
mined by the well width along with a material constant 
[19,20]. Here the components of the time-dependent po- 
sition operator of an electron in a quantum well read 
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In much the same way as above, these expressions lead to 
a zitterbewegung of electronic wave packets in a direction 
perpendicular to their group velocity. If both Rashba and 
the Dresselhaus spin-orbit coupling are present, the direc- 
tional dependence of the zitterbewegung is more compli- 
cated, which can be understood terms of the anisotropic 
dispersion relations arising in this case [21]. The corre- 
sponding expressions are rather lengthy and shall not be 
detailed here. However, in the case where the Rashba co- 
efficient is tuned to be equal in magnitude to the Dressel- 
haus term, a = ±/3, the zitterbewegung is absent. This 
is due to the additional conserved quantity which arises 
at this point and cancels the many effects of spin-orbit 
coupling [22,21]. 

Finally we mention that similar expressions can be de- 
rived for the case of heavy holes in the p-type valence 
band of III-V semiconductors (as opposed to s-type elec- 
trons studied so far) being subject to spin-orbit coupling 
due to structure-inversion asymmetry [23,24]. Again, the 
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zitterbewegung of a wave packet with its spin pointing 
initially in the z-direction is perpendicular to the group 
velocity. 

In conclusion, we have studied zitterbewegung of elec- 
tronic wave packets in III-V zinc-blende semiconductor 
quantum wells in the presence of spin-orbit coupling of 
the Rashba and Dresselhaus type. Our results suggest 
the possibility of a direct experimental proof of this os- 
cillatory motion due to relativistic effects, confirming a 
long-standing theoretical prediction. Similar results can 
be derived for the case of heavy holes in quantum wells 
under the influence of Rashba spin-orbit coupling. A very 
promising route for such experiments are high-resolution 
imaging techniques developed recently [16,17]. If the spin 
of the electron is initially aligned along the ^-direction, 
the zitterbewegung is always perpendicular to the group 
velocity of the wave packet. For possible experiments 
quantum wires are particularly attractive. For this case 
we find a resonance condition maximizing the zitterbe- 
wegung. This resonance can be reached by either tuning 
the Rashba coupling or the electron velocity along the 
wire. 
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the Swiss NSF, DARPA, ARO, ONR, and the EU Spin- 
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FIG. 1. Zitterbewegung of an electron in a harmonic quan- 
tum wire perpendicular to the wire direction. The wave num- 
ber ko y for the motion along the wire is ko y X = 5. The am- 
plitude of the zitterbewegung is maximal at the resonance 
2ako y = Tbu) (middle panel). 
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FIG. 2. Amplitude of the zitterbewegung perpendicular to 
the wire direction as a function of fi/w = 2ako y /tkj for dif- 
ferent values of the wave number ko y along the wire. 



5 



